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Capillary attraction between identical millimeter sized spheres floating at a liquid-air interface 
and the resulting aggregation is investigated at low Reynolds number. We show that the measured 
capillary forces between two spheres as a function of distance can be described by expressions 
obtained using the Nicolson approximation at low Bond numbers for far greater particle sizes than 
previously assumed. We find that viscous hydrodynamics interactions between the spheres needs 
to be included to describe the dynamics close to contact. We then consider the aggregates formed 
when a third sphere is added after the initial two spheres are already in contact. In this case, we find 



that linear superposition of capillary forces describes the observed approach qualitatively but not 
T-H quantitatively. Further, we observe an angular dependence of the structure due to a rapid decrease 

of capillary force with distance of separation which has a tendency to align the particles before 

contact. When the three particles come in contact, they may preserve their shape or rearrange to 
^ I form an equilateral triangle cluster - the lowest energy state - depending on the competition between 

^ I attraction between particles and friction. Using these observations, we demonstrate that a linear 

particle chain can be built from frictional particles with capillary attraction. 

G^ PACS numbers: 47.55.nb, 68.03.Cd, 47.55. Kf 
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I. INTRODUCTION 

Aggregates can be observed to form in particulate sys- 
tems with attractive interactions with shapes that de- 
pend on the nature of forces between particles and on 
their initial positions. An important example is the ag- 
gregation of floating objects at the surface of a liquid 
due to capillarity. This phenomenon can be easily ob- 
served in everyday examples such as clustering of bub- 
bles in a sink [T] , clumping of breakfast cereals floating 
in a bowl filled with milk [2], and biomaterials such as 
pollen or eggs of some insects species observed floating 
at the surface of ponds [3], and even swimming nema- 
todes [3] . Capillary aggregation has many important ap- 
plications as in flotation processes in ore extraction, and 
self assembly of micron-sized floating particles to fabri- 
cate new 2D-materials [5]. Further, this phenomenon has 
been exploited to study formation of ramified fractal ag- 
gregates [6], and recently by our group to examine the 
heterogeneous nature of cohesive granular media using 
spheres floating at a liquid-air interface |3- While float- 
ing particles of the same kind always attract each other, 
particles with different wetting properties can repeal each 
other. For simplicity, we limit our discussion in the fol- 
lowing to identical spherical particles. 

The first attempt to introduce capillary forces due to 
a liquid interface between floating bodies, was given by 
Poynting and Thomson 8J, who derived the forces be- 
tween semi-immersed plates. The level of fluid in between 
is increased if the plates are hydrophilic or decreased if 
the plates are hydrophobic. In both cases the curvature 
of the interface modifies hydrostatic pressure between 
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the plates, which overcomes the pressure on the exter- 
nal sides, leading to an attraction between the plates. 
This mechanism cannot be directly applied to floating 
particles [5], because gravitation, buoyancy and capillary 
forces set the height of the objects. Because the location 
of the contact line and the shape of the meniscus results 
from a competition between gravity and capillarity, its 
horizontal extension is of order of the capillary length 
= \/ 1 /{P9)i where 7 is the surface tension between 
the liquid and the atmosphere, p the density of the liquid, 
and g the gravitational acceleration. Thus, a particle in 
range of the meniscus caused by other particles, is out of 
equilibrium because the horizontal projection of surface 
tension and hydrostatic pressure is no longer isotropic. 
However, the exact calculation of capillary interactions 
becomes difficult even for spherical particles due to ge- 
ometrical complexity. Therefore, an analytical approach 
becomes possible only by simplifying the problem. 

Nicolson IJ first proposed a linear superposition ap- 
proximation to calculate the force between two identical 
bubbles at a liquid interface, which has since been ap- 
plied to fioating spheres [2 Using this approxima- 
tion implies that the calculations are restricted to small 
deformations of the interface and to particle of radius 
R which is small compared to Lc- The Bond number 
Bo = Rj^ / is the corresponding dimensionless param- 
eter to compare particle size and the capillary length and 
consequently linear calculations are limited to Bq <C 1. 
The force also depends on the contact angle 9 at the con- 
tact line between the atmosphere, the fiuid and the par- 
ticles. A complementary approach starting with Young- 
Laplace equation in bipolar coordinates for small defor- 
mations and particle sizes has been also used to derive 
the shape of the liquid interface and the force of attrac- 
tion [T0HT2] . Singh et al. [3l[T3j have developed a numeric 
simulation to study the motion of floating bodies coupled 
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with free surface flow, but the results were not corapared 
with experimental measurements. Capillary forces be- 
tween partially submerged spheres has been experimen- 
tally investigated but the particles were constrained from 
moving [14,, 15, , or with sub- millimeter sized particles at 
liquid-liquid interface where Nicolson approximation is 
expected to work JW\. 

For particles larger than a millimeter, the approxima- 
tions used in theoretical derivations discussed above be- 
come less obvious. The capillary force between millime- 
ter sized spheres has been measured by Camoin et al. [17] 
and a decreasing exponential shape was found. But the 
particles in that study were not free to move and no com- 
parison with theory was presented. A further important 
factor which has received little attention during the ag- 
gregation is the surface friction of the particles which can 
become important when more than two particles aggre- 
gate Therefore, careful experiments are necessary to 
clarify the physics of capillary aggregation for millimeter 
sized particles. 

In this paper, the force of attraction between floating 
spheres is investigated by measuring and analyzing tra- 
jectories of identical spheres with friction. In particular, 
we investigate if the mechanism of aggregation that have 
been calculated and tested for i? << Lc can be extended 
to i? ~ Lc- We will also examine dynamics of three parti- 
cles to study the effect of friction and validity of linear su- 
perposition of capillary forces for many particle systems. 
We introduce first the theoretical background related to 
the capillary interactions and the effect of viscous drag in 
Sec. m and then describe the experimental apparatus in 
Sec. 
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The experimental study of the dynamic of two 
initially isolated particles and the corresponding analysis 
are presented in Sec. |IV[ experiments with three particles 
are discussed in Sec. |V[ We conclude by demonstrating 
fabrication of a linear chain of spheres by exploiting fric- 
tion, and some remarks on the general implication of our 
study. 



II. BACKGROUND 

To discuss the nature of approximations, we first de- 
scribe in brief the derivation of the force of attraction 
between two identical floating spheres using the Nicolson 
approximation [U El |9] . Then, we consider the viscous in- 
teractions that need to be calculated to describe particle 
dynamics near contact. 



A. Capillary attraction between floating spheres 

A schematic of spheres with density Ps floating at a 
liquid-air interface and its deformation is shown in Fig. [T] 
We denote the height of the contact line relatively to the 
fluid level by Zc, and the interfacial slope at the contact 
line by z'^. The vertical position of the sphere is param- 
eterized by the angle (/)c between the lowest point of the 




FIG. 1. A schematic diagram of two spherical particles float- 
ing at a liquid-air interface. 



sphere and the level of the contact line. From geomet- 
rical considerations this angle is related to the contact 
angle by the relation: 

(t)c — — + arctan(z^) 

To find <j)c^ we use the vertical force balance for a particle 
at rest at a liquid-air interface 



P + A + T^^Q . 



(1) 



Where, P is the sphere weight, A the buoyancy for a 
semi-immersed body, and the surface tension inte- 
grated along the contact line which by symmetry is along 
the vertical direction and are given by 



A ^ttR^ pg (1 - cos^ 
and 

= 2TTjRam( 



cos 0r sm 



'-{i + z'/y/^ ■ 

Assuming that ^ 1 and Zc <C R, the following expres- 
sion can be obtained [H [5] : 



sinchc z' = sin0c tani d>c + I 



(2) 



Here, E = ^^^^ cos 6+ - cos^ 9. Given the phys- 

ical parameters of the particle and fluid, the value of 0c 
and z'^ = tan(0c + can be then numerically computed. 
Because of axisymmetry, z as a function of the distance 
to the particle center r is given by the Young-Laplace 



relation. 



z{r) 

Lr 



C{r), where C{r) is the local interface 



curvature. For small z'(r), we then obtain 
z{r) 



V^z(r) 



(3) 



This equation with boundaries conditions z — > as 
r (X) and z'(r — i?sin(0c)) — admits solution with 
modified Bessel function of the first kind Kn (x) , and fur- 
dKoix) 



ther noting that 



dx 



z{r) = -ta.n{(f)c + 9) L 



Ki{x), one obtains 



(4) 
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The numerical value of the deformation at a point along 
the contact line (meniscus size) is found by replacing 
the value of r by i? sin 

Next, we use the Nicolson approximation [T.'S^ to com- 
pute the force experienced by a particle due to the menis- 
cus of a second particle (see Fig.[T]). In this case, the hor- 
izontal projection of the capillary force integrated along 
the contact line has a resultant which is directed towards 
the center of the second sphere. The expression remains 
difficult to solve. However, using linear superposition, 
and noticing that for small deformations, the resultant 
of surface tension is modified in direction but not in am- 
plitude, Fcap can be approximated as: F^a-p = T z'2, where 

, = is the slope of interface created by the second 

^ dl r-i 
particle [18]. Using Eq. (fTj) and expressing the elemen- 
tary work of the capillary force, we obtain: 

SW = ~F,apdl={P + A) dz2. 

By integration and using Eq. (H, the energy of in- 
teraction, E{1) between these particles is expressed as 
the product of an effective weight Pg// — {P + A) — 
27r 7 i?i?oE and the liquid surface deformation z{r — I). 
Using the expression of S given by Eq. ^ , the effective 
weight can be expressed as: 

P^ff = 2n^R sin($c) tan(<I>c + 8) 



Then the energy of interaction E{1) 
written as: 



Peff z{r = I) is 



E{1) = -27r7i? 



sin($,)tan($e + g)"i^ci^o(Vi^c) 
Ki{R sin($,)/L,) 



Substituting in Fc{l) = —dE/dl, and using Eq. Q, 
we obtain the following expression for the capillary force 
between two particles: 



FS)^-CvmK^{1/L,) , 



(5) 



, , (tan(0c + 6'))2 
where, Cvm = 2n j Rsm{4)c)-rj^-—--—— . For 

Ki[R sm(0c)/^c) 
I ^ Lc, we can use the asymptotic form of Ki : Ki{x) sa 

for X 3> 1. Thus Fc decreases rapidly for dis- 
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tances larger than the capillary length. F^ is also attrac- 
tive and this feature can be explained following Singh 
and Joseph [inj. For light hydrophilic particles, rise of 
liquid between particles decreases the slope of liquid in- 
terface at the contact line. The horizontal projection of 
the tension force is therefore increased between particles, 
exceeding those on external side and leading to an attrac- 
tive force. The argument remains the same for heavier 
hydrophobic particles, by decreasing the slope of liquid 
interface at the contact line. 

It may be noted that a similar expression valid for 
small Bo is found using a different method [T0Hl2l [T6] , 
but with a slightly different prefactor: 



^ . r (tan(0c + g))^ 
with, Cpa = 27r7Lc . . 

(Ai(i? sin{(j)c)/Lc)y 

For R <^ Lc, the prefactors can be simplified, and one 

can replace the Bessel function by its asymptotic from 

Ki{x) w 1/x for a; < 1. Both Eqs. ^ and ^ then 

simplify to: 



F,o = -CoKi(l/L,) , 



(7) 



with Co = 2 7r 7 i?i3o^/^ S^. It is important to note that 
in this equation, the factor before the Bessel function was 
mainly determined by the hypothesis <C 1, whereas 
the condition of validity of Eq. ^ is given mainly by 
<C 1. Moreover in using the Nicholson approxima- 
tion, contribution of hydrostatic pressure is completely 
neglected which may not be true for millimeter sized 
spheres. While the work of AUain and Cloitre for cylin- 
ders [18] shows that pressure contribution becomes neg- 
ligible when Bo < 10, an estimate of its amplitude has 
not been reported for millimetric spheres when Bo is not 
small. Finally, it can be noted that the expression for the 
capillary force between two spheres expressed in Eqs. ([5| , 
Q and ([7]) is a product of a Bessel function —Ki[l/Lc) 
giving the spatial dependency and a constant depending 
on the model [Cvm, Cva and Cq). In the following the 
constant is simply labelled C, regardless of the theoreti- 
cal model. 

These results can be also applied to the case with more 
than two particles. Because equation Eq. ([s]) is linear, de- 
formation felt by one particle is the sum of those created 
by the other particles individually. The resulting capil- 
lary potential energy is again obtained by the Nicolson 
approximation, i.e. by multiplying the effective weight 
by the fluid surface deformation. Consequently the cap- 
illary force felt by one particle is also the sum of the force 
created by the other particles and computed using Eq. [5] 



B. Hydrodynamics of floating spheres 

Because a floating sphere begins to move due to cap- 
illary interaction, one has to consider additional hydro- 
dynamic interactions. We limit our analysis to the case 
of small Reynolds number Rf. = to neglect inertial 



to ne- 



effects and also small Capillary numbers Ca 
gleet motion of the contact line on the particle, where 
i is a typical length scale, V a typical velocity, and 
the dynamic viscosity. The drag force for a partially im- 
mersed sphere is given by the Stokes law corrected by a 
drag coefficient fd ■ 



-G-KfiRfdV, 



(8) 



F,,{l)^-CpaKi{l/L,) 



(6) 



where, v is the particle velocity, fd depends on the ver- 
tical position of the sphere relative to the interface, the 
contact angle, surface tension and density of particle and 
liquid. In order to take into account the hydrodynamic 
interactions due to the flow created by the second sphere. 
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we adopt the concept of hydrodynamic mobility intro- 
duced by Batchelor [19] for colloidal particle motion. For 
low Re, the difference of the capillary and drag forces 
projected along the axes between particles centers leads 
to: 



c,l 



c,2 



d,2 



^2CKi{l/L,)-6TrfiRfd 



dl 
di 



dl 
dt 



CK,il/L,) 



(9) 



(10) 



This last expression is corrected by multiplying it by the 
hydrodynamic mobility G{x) in terms oi x = l/Rioi two 
spheres along the line join their centers [TO], and is given 

by muni [2D]: 



G{x) 



1 



3 

2a; 



15 

4^4 



4.46 
1000 



Therefore, the equation of the motion of a sphere is given 
by 



dl 
di 



^KgG{x)K^{xR/L,) , 



(12) 



2C 



where, Kq — . This expression predicts no hy- 

e-KfifdR 

drodynamic interactions when the spheres are far from 

each other {G{x — > oo) — 1) and ^ ~ when the spheres 

come in contact (G(2) — 0), because of lubrication cre- 
ated by the liquid between the spheres. 

These results show that important features of dynam- 
ics of capillary aggregation can be extracted from the 
analysis of particles trajectories. In section [TV] we will 
study examples with two floating spheres with Bf, ^ 1 
and compare experimental results with Eq. (12 1. 



III. EXPERIMENTAL SETUP 

Experiments are conducted in a container with an alu- 
minum frame with dimensions measuring 20.3 cm long, 
20.3 cm wide, and 2.54 cm deep, and a Plexiglas bottom 
(Fig. [2| . A glass lid is placed above to avoid evaporation 
and prevent dust from falling in. To avoid boundary ef- 
fects, meniscus is pinned at a ledge along the sidewalk 
of the container in order to minimize the meniscus and 
ideally obtain a flat surface. The container is kept hor- 
izontal to within 0.1 degrees, and the level of liquid can 
be adjusted with a syringe. The particles are observed 
from above with a CMOS 1280 x 1024 pixels camera, with 
a telecentric lens, and back lit through the transparent 
bottom boundary. Images are recorded with a frame rate 
of one image per second, and particles are tracked with 
standard algorithms implemented in Image J. Relative 
error in finding particles center is less than 2% of the 
diameter. 
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FIG. 2. A crossection of the experimental apparatus with two 
particles floating at the liquid-air interface. The liquid is filled 
up to a ledge in the sidewall to minimize particle interactions 
with the container. 



In order to obtain a sufficiently viscous and dense liq- 
uid, we use a mixture of glycerol (90% by weight) and 
distilled water. The physical properties of this liquid 
are summarized in Table |l] using the data from Ref . [3T] . 
Polyethylene (High-density polyethylene) or nylon (Ny- 
lon 6) spheres are used in our study summarized in the 
Table [nl These materials are hydrophilic and lighter than 
the liquid, and were chosen in order to form stable ag- 
gregates. The contact angle 9 needed to compare with 
calculations is estimated directly by imaging a sphere 
floating on the liquid from the side. The particles are 
first immersed deep into the liquid and then allowed to 
reach their vertical equilibrium. The contact angle is 
measured from the image, as shown on Fig. [3] and val- 
ues are reported in Table [llj Using this data, we esti- 
mate the slope of the liquid interface z'^ and the am- 
plitude of meniscus Zc, using formulae in Eqs. [3] and [2] 
Even for millimetric sized particles where Bond num- 
ber is of order one, approximations used in Section 
i.e. Zc/R <C 1 and z'^ <C 1, are satisfied. Furthermore, 
we notice that the vertical position of the sphere center 
Zcenter — Zf. + Rcos(j)c is negative, so particles centers 
are below the surface, which allows us to assume that 
the drag coefficient for semi- immersed sphere /d ~ 1. 
Also when particles are touching each other, the con- 
tact point between particles is located below the liquid 
surface, which prevents high interface distortion when 
particles are close. 



Finally, we note that experiments are very sensitive 
to vibrations and thermal convection. Therefore, room 
temperature is regulated at 23° C and several hours are 
needed after filling the container to reach thermal equi- 
librium to commence experiments. 
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Density p (kgm 


Surface tension 7 (N m 


^) Viscosity /i (Pas) Capillary length Lc (mm) 


1233 ± 3 


0.064 ± 0.002 


0.175 ±0.015 2.30 ± 0.04 



TABLE I. Physical properties of glycerol- water mixture used in the experiments [21]. 



Material 


D (mm) p 


s (kg m-^) 


6 (degrees) 


Bo 


(f)c (degrees) 


Zc (mm) 


4 


Zcenter (mm) 


Polyethylene 


6.35 ± 0.05 


950 


40 ± 3 


1.91 ± 0.1 


124 


0.49 


- 0.29 


-1.27 


Polyethylene 


3.175 ± 0.05 


950 


20 ± 2 


0.48 ± 0.03 


151 


0.16 


- 0.15 


-1.23 


Nylon 


3.175 ± 0.03 


1150 


5 ± 5 


0.48 ± 0.03 


169 


0.074 


- 0.11 


-1.48 


Nylon 


2.38125 ± 0.05 


1150 


5 ± 5 


0.27 ± 0.02 


171 


0.038 


- 0.075 


-1.14 


Nylon 


1.5875 ± 0.05 


1150 


5 ± 5 


0.12 ± 0.009 


173 


0.014 


- 0.042 


-0.77 



TABLE II. The diameter D and the density ps of the particles were given by the provider, the contact angle 8 is measured 
directly (see Fig. [3|, and the Bond number Bo is the ratio R^/Lc^. The other parameters characterize the equilibrium position 
of a floating sphere including the angle (j)c computed from Eq. Amplitude of the meniscus around the particles Zc and 
meniscus slope z'^ are estimated using Eqs. ([3| and H. Vertical position of sphere center is given by Zcenter = Zc + Rcoscpc- 




FIG. 3. A side view of a floating particle with D = 3.175 mm. 
The contact angle corresponds to the angle between the 
liquid-air interface and the tangent to the particle at the same 
point. 



IV. ATTRACTION BETWEEN TWO 
FLOATING PARTICLES 



In this section, we discuss the motion of two initially 
separated floating particles. To obtain reproducible data, 
the particles are first totally immersed in the liquid and 
then placed in their initial position with negligible veloc- 
ities. We plot the distance between the centers of two 
polyethylene spheres I as a function of time in Fig. |4] 
Curves for various initial positions overlap with each 
other when we chose the time of origin to be the time 
of contact. This curve shows that the particle dynamics 
before contact, is a function only of separation distance. 

We then determine the velocity of the spheres w^ei = 
dl/dt from the separation distance over time and plot it 
in Fig. [5] versus x = l/R. The data was averaged over 
ten experiments to reduce statistical and measurement 
errors. The velocity is observed to first increase rapidly 
as the spheres approach each other, and then decrease as 
the particles approach each other because of lubrication, 
and reaches zero once particles come in contact due to 
hard core repulsion. 

From this measured velocity, we can extract the drag 
force experienced by the particles, and assuming that in- 
ertia is negligible, we can also obtain the capillary force 
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FIG. 4. The distance between the centers of the two polyethy- 
lene particles versus time {D = 3.175mm). The time is cho- 
sen to be equal to zero when the spheres come into contact. 
The different symbols correspond to various initial separation 
distance. 



by equating it with the drag force. Now the maximum 



velocity is approximately Vrei /2 ^ 
fore, the Reynolds number. Re = 

and the Capillary number, Ca — 



4x10" 
VRp 

A* 

VRn 

PI 



ms ^. There- 

: 4 X 10"-* < 1 

9 X 10-* < 1 



which are consistent with our assumptions. 

We fit the experimental data with Eq. (12 1, and find 



that its form is well described with Kq = 9.10 x 10~ ± 
0.1 X lO^"* ms^^ (see Fig.jsJ. If we assume that the drag 
coefficient for semi- immersed particle /d = 1, then from 

n C 

Eq. (5) we obtain Kvm = 7, ^ = 9.66 x lO^'^ms^^. 

37r/iit 

Considering the significant number of assumptions, the 
difference of about 6 % between the fit and the theory 
is remarkably small. Further, the decrease of velocity at 
small particles separation {l/R < 2.5 in Fig.js]) is well de- 
scribed by the Batchelor model of viscous hydrodynamics 
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FIG. 5. Approach velocity Vrei versus x = l/R of Polyethylene 
spheres with D = 3.175 mm. Errors bars are estimated from 
the statistical dispers ion using 12 data sets. The curve is 



fitted by equation Eq. ( 12 \ with coefRcient Kg = 9.10x 10 ± 
0.1 X 10"*m.s"\ 



interactions at low Re. 

We measured Kq and Kvm with spheres with other 
sizes and material properties listed Table |III[ The mea- 
sured velocity is scaled with Kq and compared with 
Eq. ^ inFig.|6] Good agreement is observed except for 
the largest separation for Nylon spheres. The difference 
may arise due to various reasons. First, Eq. ([t]) was es- 
tablished assuming that the deformation of the interface 
and the slope is small. Results for Polyethylene spheres 
of diameter D — 6.35 mm are less in agreement with the 
model, because we have z^/R = 0.154 and z'^ = —0.29, 
which are not so small compared with 1. Second, we 
assumed fd = 1, which may not be accurate. Finally, 
Nylon particles have a higher density, and the resulting 
capillary force is of an order of magnitude smaller and 
difScult to measure for large l/Lc- 

Overall, our experiments show that capillary attrac- 
tion between millimetric floating spheres is reasonably 
well described by the calculations presented in Sec. [ll] 
For completeness, we present the theoretical results ob- 
tained with Eq. ([6| in Tab. Ill The agreement is similar 
to those obtained with Eq. (5fexcept for the highest Kq- 



While it is not possible to be conclusive as to which model 
is more accurate, it appears that the Nicolson approxi- 
mation [211^ is simpler to apply to experiments. Further, 
our experiments show that both models remain valid for 
Bp 1 at least when the liquid surface is not strongly 
distorted even during contact as in our case where parti- 
cle centers are well below the liquid surface. 

Our results lead to interesting information on the ag- 



gregation process. First, we can estimate the aggregation 
time using the equation of motion when we neglect hy- 
drodynamics interactions. This approximation is reason- 
able because the time scale of the approach is dominated 
by the time when particles are far from each other when 
viscous hydrodynamic interaction between particles can 
be neglected (see Fig. [4]) . If one has a system of isolated 
spheres, and take the initial separation equal to the mean 
separation of the particles Irm we estimate the contact 
time: 



tn = 



2R 



-1 



■dl. 



(13) 



These experiments were performed with 10 particles 
{D = 3.175 mm), with an initial separation Im = 1.85 cm. 
We find that tc — 13000 s, where as the estimate gives 
tcilm) = 24000s which correctly captures its order of 
magnitude. 

Moreover the cohesive force inside an aggregate can be 
estimated using Eq. ([5|, with I = 2R. Because the con- 
tact point between the spheres is below the surface, there 
is indeed no discontinuity of the interface. For spheres 
with D = 3.175 mm, and using the corresponding value 
of Kg, we find the cohesive force at contact: 

Fc(2i?) = Zt: ^iRKgKi{2R/L^) w 7.9 x 10~^N, 

which is significantly smaller than the effective weight of 
a particle « 4.3 x 10^^ N. Finally the good agreement 
between measurements and theoretical results obtained 
with Nicolson approximation implies that the hydrostatic 
pressure contribution is negligible which is not obvious 
when Bo ^ 1- Force due to hydrostatic pressure re- 
sults from the difference of liquid level around the parti- 
cle caused by the meniscus created by a second particle. 
Integration of the pressure on a sphere with a spatially 
varying contact line cannot be accomplished analytically. 
Using linear superposition of surface deformations and 
expression of the capillary force between two plates [H [5] , 
an estimate of the force felt by the first particle due to 
hydrostatic pressure can be obtained as: 



Fp{2R) « pgR{z2{mY - z^iR) 



where Z2 (r) , the deformation of the liquid surface due to 
the second particle only is computed using Eq. ^. For 
D = 3.175mm, we find Fp{2R) « 1.3 x lO-'^N. For the 
largest spheres with D = 6.35 mm, Fp{2R) « 4.9xlO~^N 
when the capillary force becomes Fc{2R) w 5.6 x lO^^N. 
Consequently it appears that for particles of similar or 
larger size, hydrostatic pressure contribution should be 
taken in account to compute amplitude of capillary at- 
traction. 



V. EXPERIMENTS WITH THREE PARTICLES 

We now discuss the formation of aggregates when a 
third sphere is added at various positions relative to the 




FIG. 6. Average dimensionless approach velocity Vrei/Kc versus l/Lc, for five different kinds of particles listed in Table III] 



Measurements are compared with the spatial dependency given by Eq. ( 12 1 



Material 



D (mm) Experimental coefficient Theoretical coefficient Theoretical coefficient 







Kg (ms 






KvM (ms 






KpA (ms" 




Polyethylene 


6.35 


2.0 X 10" 


±0.3 X 10 


-2 


3.7 X 10" 


-2 


±0.7 X 10" 


-2 


6.8 X 10"'' 


± 1.2 X 10"^ 


Polyethylene 


3.175 


9.1 X 10" 


"''±0.1 X 10 


-4 


9.7 X 10" 


4 


±0.7 X 10" 


-4 


11 X 10"* 


±0.7 X lO""* 


Nylon 


3.175 


15 X 10" 


5 ±0.1 X 10" 


5 


7.9 X 10" 


5 


± 1.2 X 10" 


-5 


8.1 X 10"^ 


± 1.4 X 10"'' 


Nylon 


2.38125 


4.2 X 10" 


±0.1 X 10 


-5 


1.8 X 10" 


5 


±0.2 X 10" 


-5 


1.9 X lO"'^ 


±0.2 X 10"^ 


Nylon 


1.5875 


4.5 X 10" 


±0.2 X 10 


-6 


2.4 X 10" 


-6 


±0.4 X 10" 


-6 


2.4 X lO"'' 


±0.4 X 10"'^ 



TABLE III. The coefficients Kg observed in the experiments which sets the amplitude of capillary attraction. The experimental 
values were obtained using Eq. \12\ and errors are estimated from accuracy of the fit. The theoretical values were calculated 
from the Eqs. Q, ([5| and ([7|, the errors arise from the uncertainties in the parameters values. 



two sphere cluster. "We label particle 1 and 2 as the ones 
that are initially in contact, with the first being closer to 
the third particle. We label Im as the distance between 
the third particle, and M, the point of contact between 
particle 1 and 2, and lij is the distance between particles 
i and j. Inside the cluster, we have I12 = 2R — D. (3 
is the angle between particle 1 and 2 and particle 3 and 
M (see Fig. [?]). After the particles are placed in their 
initial positions, we observe that the cluster of two par- 
ticles moves as a rigid body and can rotate depending 
on the initial value of /3. Once contact occurs between 
the cluster of two particles and the third particle, a re- 
arrangement can occur depending on initial conditions 
producing an equilateral triangle (Fig. [s]). 



A. Capillary force due to a cluster of two particles 
on a third particle before contact 

We define Fc as the force exerted on the particle 
3 by the two particle cluster (see Fig. [t]). We have 
Fi ^ CKiihs/Lc) ui and Fa = C Ki{l23/Lc) U2, corre- 
sponding to forces which would be exerted by particles 
1 and 2 respectively on particle 3, ignoring the effect of 
the other particle, and Ui and U2 are associated unitary 

vectors. The approach velocity is defined — p-. We con- 

dt 

sider the projection of Fc along this direction using an 
unitary vector u. Then assuming linear superposition, 
we obtain: 

Fc u = Fi u + F2 • u . 

Drag forces are difficult to estimate in presence of a clus- 
ter of two particles. As hydrodynamic interactions are 
dominated by the contribution of the closest particle. 
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^ Particle 1 



FIG. 7. The angle /3 and a corresponding to the three particle 
system viewed from above. 




(a) 



(b) 



(c) 



FIG. 8. A sequence of images showing aggregation of three 
particles /3i = 77 degrees and /3/ = 90 degrees. 



hydrodynamic mobility can be estimated as G{li3/R) in 
first approximation. Then, 

- dhi/dt = Kg Gihs/R) K, (h^/L,) Ui • u 

+ Kg Gihs/R) Kiihs/L,) U2 • u (14) 

/D 1, then we have Ki{Im/Lc) ~ 

'.M/La_ Therefore, Fi and F2 decrease rapidly 



If 



'A/ 



21 

with separation distance between particles. As a result, 
the closer particle is attracted more strongly which causes 
the cluster to rotate and decreases /3. 



Next, we test if linear superposition given by Eq. (14) 
can be used to describe the formation of the cluster. The 
approach velocity Vrei is plotted versus the distance be- 
tween particle 1 and 3 Z13 in Fig. |9] for (3{t — 0) — f]^ — 
degrees (a) and I3{t = 0) — f3i — 90 degrees (b). The 
curves are quite similar to the one obtained previously 
for attraction between two particles (Fig. [5| . Then the 
right part of Eq. ( 14 ) and the part exerted by particle 
1 only, are added in Fig. [9j We observe that the force 
exerted by particle 1 dominates as Z13 < I23 but is not 
sufficient to describe the total force (Fig.lo]). When par- 
ticles are initially aligned {(3i = degrees (a)), the linear 
superposition is a good approximation. Decrease of Vrei 
close to contact due to lubrication is also well reproduced 
using mobility for the closest particle. For a larger initial 
angle Pi , result of linear superposition fails to follow ex- 
perimental curve when the third particle approaches the 
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FIG. 9. Test of the linear superposition assumption for 
/? = degrees (a) and = 90 degrees (b). Average approach 
velocity of the third particle is plotted as a function of the 
distance I. 



cluster. Failure of linear superposition may not be sur- 
prising because the deformation for two particles is not 
equal to the sum of the deformation of two individual 
spheres when the spheres are close to each other. Nev- 
ertheless, the linear superposition gives the correct qual- 
itative description, and even gives a reasonable estimate 
of Fc when a particle is at a sufficiently large distance 
from the two-sphere cluster. 
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FIG. 10. (a) Evolution of the angle /3 for polyethylene spheres 
(D = 3.175 mm) with two different initial values, 77° and 37°. 
The cluster always turns before contact. For large values of 
j3i, it also turns after contact, and this evolution is faster 
than before contact, (b) The angle between the third particle 
and the center of the cluster of two particles at the instant 
of contact /3c (red circles) and after any rearrangement as 
occurred j3f (blue crosses) as a function of its initial value Pi. 



B. Particle Friction: Shape and dynamics after 
contact 



We plot the evolution of /3 for two initial angles f3i in 
Fig. [lOl^a) to illustrate the two different kinds of clusters 
observed. In both cases, /3 decreases before contact, but 
while the newly formed 3-particle cluster stops evolving 



Material 


D 


/3c,t 




a 






k 




(mm) 


(degrees 


) (deg 


Tees) 




Polyethylene 


6.35 


41 ± 


4 


119 


± 


4 


0.0462 


Polyethylene 


3.175 


33 ± 


1 


131 


± 


1 


0.0761 


Nylon 


3.175 


21 ± 


1 


149 


± 


1 


0.0417 


Nylon 


2.38125 


29 ± 


3 


136 


± 


3 


0.0819 


Nylon 


1.5875 


35 ± 


2 


128 


± 


2 


0.124 



TABLE IV. Transition angles at contact when triangular com- 
pact cluster form, and the effective coefficient of friction for 
the various kinds of spheres investigated. 



for small initial angles, a rapid rearrangement is observed 
where /3 increases to 90 degrees for the larger initial an- 
gle. This second situation corresponds to formation of an 
equilateral triangle which is the lowest energy state for a 
floating three particle cluster. In order to investigate the 
dependence of the two final states on initial conditions, 
experiments were performed systematically as a function 
of f3i and are represented in Fig. 10 'b) by plotting the 
angle /3c at contact and /3j after any rearrangement has 
occurred. 

As argued previously, the rapid decrease of force with 
distance tends to decrease f3 before contact as the spheres 
approach each other. For the smallest and largest values 
of /3i, cluster do not rotate significantly and /3c remains 
very close to /3i. Then after contact, if /3i is smaller than 
Pi^t = 62 ± 2 degree, the shape of the three particle 
cluster does not evolve, and thus /3/ = /3c. Otherwise 
capillary attraction inside the cluster modifies its shape 
and tends to increase the final angle to 90 degree. This 
shape corresponds to a regular triangular cluster where 
all the spheres are in contact with each other. As the 
distance between particles are minimal, this state has 
the lowest possible energy. 

We explain the difference of behavior because of the 
presence of friction which prevents the spheres from 
rolling and sliding on each other once they come in con- 
tact. Figure 10 h) can be considered as a transition dia- 
gram between two phases: compact and aligned clusters 
due to the competition between friction and capillarity. 
Below a certain value I3i^t of /3, final cluster shape is 
aligned, and the final shape is compact above this value. 
The few exceptions observed may be due to particle sur- 
face irregularities like small bubbles or dust trapped on 
the spheres. 

It appears pertinent to express the transition at the 
contact instant, which gives f3c,t = 33 ± 1 degree and 
relate this value to the angle a between I12 and ^13 
(see Fig. [t]). The transition occurs when a is below 
Q!c,t — 131 ± 1 degrees, when friction can no longer bal- 
ance the capillary force. Similar behavior is observed us- 
ing other kinds of particles with different transitions an- 



gles ac,t which are reported in Table IV Using the tran- 
sition angles etc,*, one can estimate the effective friction 
coefhcicnt k between particles using the ratio of tangen- 
tial and normal component of force acting on particle 3 at 
the point of contact with particle 1. The tangential com- 
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ponent is entirely due to force ^23 due to particle 2, and 
is given by F23 cos(ac,t/2). The normal component is due 
to the sum of force F13, and component F23 sin(Q;c,t/2) 
parallel with ^"13. Therefore, 

^ _ F23 cos(ae,t/2) 

F13 + F23 sin(ac,t/2)' 

Because F13 = CKi{D/L^) and F23 = CKiihs/L^), we 
obtain 

, ^ Ki{l23/Lc)cos{ac,t/2) 
- K, {D/L,) + K, {I23/L,) sin(a,,t/2) ' ^ ^ 

Where, I23 = D^/2{1 — cos Qfc.*) when particles 1 and 3 
are in contact. Computed values of k are reported in Ta- 
ble |IV| It may be noted that k decreases with particle size 
and further depends on the material. While rolling and 
sliding friction coefficients for these materials on various 
hard surfaces have been reported and depend on normal 
force [? ], we were unable to find corresponding data to 
compare to when particles are immersed in a liquid which 
may be expected to modify the friction at contact. 

C. Linear chain assembly with capillary and 
friction forces 

Building on our observation, we demonstrate an exam- 
ple of an aggregate with frictional particles self-organized 
with capillary forces which is not possible with friction- 
less particles. The chain is grown initially one particle 
at a time in Fig. |ll| [a-d) by dropping a particle near one 
end of a the developing chain. Then, the chain can be 
straightened to remove any kinks by introducing particles 
near opposite ends which tends to straighten the chain 
(see Fig. [nje)). The images shown here correspond to 
polyethylene spheres with D = 6.35 mm. Using this tech- 
nique, we were able to grow chains as long as 12 particles 
as shown in Fig.[Tl|f). Longer chains tend to be unstable 
and fold to form two rows of particles which is energeti- 
cally more favorable. 

VI. CONCLUSIONS 

We have experimentally investigated capillary aggre- 
gation of small number of floating spheres and shown 
that attraction between two spheres is reasonably well 
described by expression of capillary force originally de- 
veloped for small Bond number Bo |2]. Specifically, we 
have shown that for particles which are hydrophilic and 
lighter than the fluid, interface deformation and its slope 
remain small even when Bond number is of order one, 
allowing one to extend the range of validity of the the- 
ory. Further, our measurements are sufficiently accurate 



to also demonstrate the crucial role of hydrodynamic in- 
teractions producing rapid decrease of spheres velocity as 
they come in contact due to viscous efl^ects. 



(a) 


• 




1*^ 


(b) 


m •••• 


(e) 




(c) 









FIG. 11. (a-f) An example of a chain of twelve spheres ob- 
tained by adding successively particles along the cluster axis. 
The particles are assembled using capillary forces and stabi- 
lized by friction. 

We have also studied aggregation of a third particle 
after formation of a cluster of two particles. The theo- 
retical estimate using linear superposition appears to be 
less accurate in this case, but can be still used to ex- 
plain the observations qualitatively. The final shape of 
cluster is determined by the initial position of the third 
particle relatively to the cluster of two particles. Rota- 
tion of the cluster, due to the differential attraction of 
the two sphere tends to create aligned clusters of three 
particles. But if the angle /? between the three particles 
is small enough, capillary forces between the two parti- 
cles at the ends overcomes friction to form a compact 
triangular cluster. 

Finally, we note that the observation of large heteroge- 
neous structures in larger capillary aggregates [7j can be 
also explained by our analysis. Porous clusters were ini- 
tially observed to evolve in that study from dilute concen- 
tration of particles. As illustrated in our study, particles 
will rearrange to form a compact triangular cluster only 
if friction at contact can be overcome by capillary forces. 
Thus chain like structures which connect various parts of 
the aggregate can be formed that encompass pores which 
are locked in place unless compression is applied. Even as 
particle concentration is increased, such regions lead to 
defects in the more dense aggregates that form. Our ob- 
servations demonstrate that friction is clearly necessary 
to explain shapes of capillary aggregates at the millimeter 
scale. 
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